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A mean-field treatment of the phase string effect in the t — J model is presented. Such a theory 
is able to unite the antiferromagnetic (AF) phase at half-filling and metallic phase at finite doping 
within a single theoretical framework. We find that the low-temperature occurrence of the AF long 
range ordering (AFLRO) at half-filling and superconducting condensation in metallic phase are all 
due to Bose condensations of spinons and holons, respectively, on the top of a spin background 
described by bosonic resonating-valence-bond (RVB) pairing. The fact that both spinon and holon 
here are bosonic objects, as the result of the phase string effect, represents a crucial difference from 
the conventional slave-boson and slave-fermion approaches. This theory also allows an underdoped 
metallic regime where the Bose condensation of spinons can still exist. Even though the AFLRO 
is gone here, such a regime corresponds to a microscopic charge inhomogeneity with short-ranged 
spin ordering. We discuss some characteristic experimental consequences for those different metallic 
regimes. A perspective on broader issues based on the phase string theory is also discussed. 
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I. INTRODUCTION 



The t — J Hamiltonian is one of the simplest nontrivial model to describe how doped holes move on the AF spin 
background and is widely used to characterize the physics in the CuOi layers of cuprates. A tremendous effort 
has been contributed to the investigation of the t — J model. The most popular approaches to the metallic phase 
are often based on the so-called slave-boson method |l| in which the degrees of freedom associated with spins are 
described in terms of fermionic description. There have been many proposals of mean-field ground states based on 
such a fermionic description of spins, ranging from the earlier fermionic RVB states, Q gauge-theory description, 
HQ SU(2) formalism, Q to possible fractional-statistics. |(|f?J However, there is an inherent problem quite general to 
the fermionic description of spins: approaches based on it usually fail in faithfully producing correct AF correlations 
especially at small doping. 

At half-filling, for example, the exact ground state is known to satisfy the Marshall sign rule ^ (for a bipartite 
lattice) but a fermionic description of spins would show redundant signs: even exchanging two same spins will give 
rise to a sign change of wavefunction due to the fermionic statistics. Under strict enforcement of no double occupancy 
constraint, those unphysical signs would not have any effect. But in mean-field approximations, this "sign problem" 
will always show up and cause serious problem like overall underestimate of AF correlations. 

By contrast, the Marshall sign can be easily incorporated into a bosonic description of spin degrees of freedom, where 
no extra sign problem would be caused by the statistics of bosons. It is the reason contributing to the success of the 
bosonic RVB description PJTo|] and its mean-field version - the Schwinger-boson mean-field approach pd[ | in describing 
spin properties at half-filling. A variational wavefunction based on the bosonic RVB picture can produce p|]lO|] an 
unrivaled accurate ground-state energy (— 0.3344J per bond as compared to the exact numerical value of — 0.3346J 
per bond for the Heisenberg model), and a generalized description [fiof can precisely provide not only the ground-state 
energy, staggered magnetization, but also spin excitation spectrum in the whole Brillouin zone. Therefore, a bosonic 
description of spins seems most natural at half-filling. 

It is thus tempting for one to use the bosonic description of spins as a starting point and try to get into the metallic 
phase by doping. The connection between antiferromagnetic and metallic phases is commonly perceived important 
in the t — J model, and is also believed by many to be the key in search for the mechanism of superconductivity 
in cuprates. Unfortunately, the mean-field study in the Schwinger-boson-slave-fermion approach Jl2]of the t — J 
model, which is based on the bosonic description of spins and has been quite successful at half- filling, [|ll| soon meets 
problematic consequences once holes are introduced - encountering the so-called spiral phase and its derivatives. 
IjjIjijJ It seems that one could not avoid such a spiral instability so long as a perturbative approach is adopted. (fif l 
One of many problems with spiral phases involves an underestimated kinetic energy tx 5 2 ^j at weak-doping S <C 1, 
which is also accompanied by a very quick descent to ferromagnetic phase at slightly larger S. p4| 

This implies that doped holes may have introduced some singular doping effect which has been mistreated in the 
mean-field approximations. Such a singular doping effect has been recently identified [|lft[17| by reexamining the 
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motion of doped holes in the AF background. It has been found that spin-mismatches caused by the hopping of 
doped holes cannot be completely "repaired" through spin flips at low energy. Such a residual nonrepairable effect 
can be expressed by a path-dependent phase product known as phase string. [|l6| Due to the phase string effect, a 
hole slowly moves through a closed path will acquire a nontrivial Berry's phase. As this phase string effect is very 
singular locally at a lattice constant scale, its topological effect can be easily lost if a conventional mean-field average 
is involved — a reason causing the aforementioned spiral instability. 

In order to handle such a singular phase string effect hidden in the conventional Schwinger-boson-slave-fermion 
scheme, a unitary transformation |l7[ has been introduced to regulate the Hamiltonian such that the local singularity of 
the phase string (at the scale of one lattice constant) is "gauged away" , while its large distance topological consequence 
is explicitly incorporated into the Hamiltonian. The resulting exact reformulation JlTj of the t — J model is believed 
to be more suitable for a perturbative treatment, in contrast to the original slave- fermion formalism. The underlying 
physical implication is that the "holon" and "spinon" defined in the slave-fermion scheme |l2|,|l4| may not be really 
separable due to the hidden phase string effect, but those in the new formalism may become truely elementary 
excitations. In the one-dimensional (ID) case, correct Luttinger-liquid behaviors indeed can be reproduced |P"t|| after 
a mean-field decoupling of the spin and charge degrees of freedom in this new scheme. 

In this paper, we develop a generalized mean-field-type theory based on this new formalism of the t — J model 
in the two-dimensional (2D) case. This theory recovers the well-known Schwinger-boson mean-field state |lf| at 
half-filling while predicts a metallic phase at finite doping without encountering any spiral instability. It offers a 
unified phase diagram for the t — J model at small doping, in which an insulating AFLRO phase, an underdoped 
metallic phase with the phenomena of pseudo-gap and charge inhomogeneity, as well as a uniform metallic phase 
with "optimized" superconducting transition temperature, are all natural consequences happening on a single spin 
background controlled by bosonic RVB order pairing. The phase string effect plays a crucial role here to connect those 
different phases together within a single theoretical framework. A short version of this work was published earlier. 




II. MEAN-FIELD THEORY BASED ON PHASE STRING EFFECT 



In the standard slave-fermion formalism of the t — J Hamiltonian, electron annihilation operator Ci a is written as 

Cia = flb la (-a)\ (2.1) 

in which fj is fermionic "holon" creation operator and is bosonic (Schwinger-boson) "spinon" annihilation operator, 
satisfying no double occupancy constraint 



fi fi 



(2.2) 



The t — J model, H t -j = H t + Hj, is composed of two terms: the hopping term H t is given by 



Ht 



H.c, 



and the superexchange term is 



Hj 



(2.3) 



(2.4) 



Note that the staggered pha se factor (— a) 1 in (2.1) is introduced jl7) to explicitly tra ck the Marshall sign, which leads 
to the negative sign in ( |2.4| ). A sign a = ±1 then appears in the hopping term ( [2.3[ ) which is the origin of the phase 
string effect [ pl||ri| mentioned in the Introduction. Due to such a sign, a hole moving from a site a to an another 
site b will acquire a sequence of signs, i.e., a phase string as shown in Fig. 1, which has been shown [jl6 17 to be 
nonrepairable by the spin flip process governed by Hj . It implies that the slave-fermion formalism of the t — J model 
cannot be treated in a perturbative way in doped case. 
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A. Phase string representation 



It has been shown that the above singular effect of phase string can be regulated after a unitary transformation. 
Jl7[ The resulting new formalism is known as the phase string representation. The hopping term H t in this new 
representation becomes [ fl7| 

H t = -t (e<) h\h 3 (e-4«) b \ a b ir , + H.c. (2.5) 
and the superexchange term Hj reads 

J 

(ij)aa 



Hj = ~ E b l b U ('" T ' r .) • (2-8) 



Note that the fermionic operator fa now is replaced by a bosonic holon operator hi in this new formalism. So one role 
of the phase string effect is to turn holons from fermions into bosons. Both holon and spinon are now described by 
bosonic operators which still satisfy the following no double occupancy constraint 



In this new formalism, the singular phase string effect, as represented by the sign a in the hopping term ( |2.3| ) of 
the original slave-fermion representation, is "gauged away" , but its topological effect is left and exactly tracked by 
lattice gauge fields A{j and . These fields are defined as follows 

4=4?.-$ (2.8) 

with 



(2.9) 



4 = jE^(0-%(0], (2-10) 

and 
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4 = ^E^W-^(0]nf. (2.11) 
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Here n\ a and are spinon and holon number operators, respectively. #,(/) is defined as an angle 

9i(l) = Im In ( Zi - zi) (2.12) 

with Zi — Xi + iyi representing the complex coordinate of a lattice site i. 

The physical meaning of Afj and A^ have been discussed in Ref . [0 : Afj and A^ describe quantized flux tubes 
bound to spinons and holons, respectively (Fig. 2 illustrates the case for A^A. Furthermore, the field <p^ describes a 
uniform flux threading through the 2D plane with a strength ir per plaquette: Y] n 0? = ±7r. It is also noted that a 



7r-flux neutral topological excitation has been previously discussed 19 in the pure Heisenberg model, which resembles 
a quantized flux line in the mixed phase of a BCS superconductor. Here the flux quanta are bound to the doped holes 
due to the phase string effect. 

Electron operator in this representation becomes |L7| 



h}b la (-aYe^ . (2.13) 

mlocal phase factor e ie ^ 9 precisely 

defined by 



Here the nonlocal phase factor e 10 ^ 9 precisely keeps the track of the singular part of the phase string effect and is 
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with 



K Hn9 = \[*\-<>S% (2-14) 



$ t > = 5>(0 E . ( 2 - 15 ) 



and 



S^E^K- (2.16) 



B. Mean- field approximation 

For the sake of clarity, in the following we first consider the superexchange term Hj and then include the hopping 
term Ht, due to different natures represented by them. 

1. Generalized mean-field treatment of Hj 

At half-filling, the mean-field theory based on the bosonic RVB picture is known as the Schwinger-boson mean-field 
theory which was first introduced by Arovas and Auerbach |Tl| . Such a mean- field is characterized by a bosonic RVB 
order parameter 

A S = ^(hrbj-r) (2.17) 

(7 

for the nearest-neighbor sites i and j [i = nn(j)]. 

The present formalism only differs from the Schwinger-boson, slave-fermion formalism in doped case, where a gauge 
field Ah emerges. Since spinons are subject to this gauge-field in Hj, it is natural to incorporating the link 

variable e~ laAi i into the order parameter ( 2.1 7| ) . Namely, 



(7 

Afj defined here is then "gauge- invariant" under an "internal" gauge transformation: A'^ — > A\= + 6i — 6j, and 
b — > b- p icrSi 



Based on such an order parameter, one may write down the mean-field version of Hj in (2.6) in a standard procedure 



H J a = — E K-)* e-^hrb^ + H.c. ' ^E A -. " ' A (E b « b *° - C 1 - V N ) ■ (2 1 !) 1 

(ij)cr (ij) \ ia I 



where the last term with a Lagrangian multiplier A is introduced to enforce the condition of total spinon number, 
with TV denoting the total lattice number and S doping concentration. In order to diagonalize Hj , we introduce the 
following Bogoliubov transformation 



bier = E [Urnaiihrncy ~ *W(i)7m-a) ! (2.20) 
m 

and seek the solution 

[Hi, 1 l a ]=E ml l a , (2.21) 

[Hi^ m<T )=-E mlm<T . (2.22) 
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Here 7 mcr and 7^ are bosonic annihilation and creation operators, respectively, for an eigenstate with quantum 
number m and spin a. In terms of bosonic commutation relations, one easily finds that u ma (i) and v mcr (i) satisfy 



and 



According to ( 2.19| ), we have 



E [Uma{i)u* ma (j) - V ma (i)v* ma (j)} = S i: 



E [u m a(i)v m -a{j) - V ma {l)u m - a {j)} = 0. 



j=nn(i) 



Then by using (2.20), (2.21) and ( [2.22 ), it is straightforward to derive the following relations from (2.25) 



^m^mtr(^) 



J 



E A? je -^<v* m _ a (j) - Xu mrT (i), 



j=nn(i) 



J 



E A« ie -^<_ CT (i) + A« mCT (i). 



j=nn(i) 

We can further express u ma {i) and v ma -(i) by the "one-particle" wavefunction w ma (i) as follows 

^mer(^) — ^m^mcr(^); 

where u m and u m will be taken to be real and satisfy 

Then w ma -(i) is normalized following from ( [2.30| ): 

^w m(T (i)<<r(i) = 5 H- 

m 

Equations ( |2.26 ) and (2.27) then reduce to an eigen-equation for the one-particle wavefunction 



j=nn(i) 



The eigenvalue £ m in ( [2.32 ) is related to E m and u m and w m as follows 



{E m + A m )- 



(2.23) 
(2.24) 

(2.25) 

(2.26) 
(2.27) 



(2.28) 
(2.29) 



(2.30) 



(2.31) 



(2.32) 



(2.33) 



Here A m is the same as the Lagrangian multiplier A, but we write it in a general form because later it will be modified 
once the hopping term is introduced. In terms of (2.33) and (2.30), one obtains 



Em — V^m Cm) 



and 



V2\E m J 



1/2 



1/2 



(2.34) 



(2.35) 



(2.36) 
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The signs of u m and v m are determined up to sgn(v m /u m ) — sgn(£ m ). As a convention we will always choose 
u m = \u rn \ and v m = \v m \sgn(^ m ). 

Thus Hg is diagonalized as = J2 m(J -E m 7 m(7 7 mcr + const , according to ( 2.21 ) and ( 2.22j ). The order parameter 
Afj can be self-consistently determined by the definition ( 2.18| ) as 



A? 



E< 



-to A 



1 +X!(7ma7r7 



(2.37) 



In the following we will always consider the solution of a real order parameter Ay. In this case, it can be checked 
self-consistently that w ma = w* m _ a according to (2.32) and (Af^) = Ay in terms of (2.37). The order-parameter 



equation may be further simplified if one multiplies ( |2~37| ) by (Ay) and sums over (ij) with using (|2.32[): 



(ij) 



JE n 



(2.38) 



with (7^ CT 7 mCT ) = l/(e /3£ '" - 1) and = l/k B T. Finally, the condition 



(E 6 tM = (i-w 

which is enforced by the Lagrangian multiplier in ( 2.19| ) and can be rewritten as 

2 -<^4E 



Am 0E 
m#0 m 



(2.39) 



(2.40) 



where n b BC is introduced to describe the contribution from E rn — state (denoted by m — 0) when the Bosc 
condensation of spinons occurs. |po| In comparison with the zero-doping Schwinger-boson mean- field theory, the 
above Bogoliubov-de Gennes scheme at finite doping mainly differs in the one-particle eigenequation (2.32) (and 
the resulting energy spectrum £ m and wavefunction w ma -(i)). One may simply recover the Schwinger-boson mean- 
field results by setting Ay — in (2.32) and obtaining = — JAg(cos k x a + cosk y a) and the Bloch wavefundtion 



Who 



i e iakn m ^- ne n0 _j 10 i e case. 



So the doping effect has entered the above mean-field theory in two ways: one is through the particle number 
condition in ( 2.39| ); the other is through the gauge field Ay. For the mean- field theory to work, Ay has been 
implicitly assumed as a time-independent field. But with holons moving around, A\^ will usually gain a dynamic 
effect. To see this, let us consider the following gauge-invariant quantity 



E4=-E^ 



c tec 

where C is an arbitrary counter-clockwise closed path. If one redefines 

n\ = 5 + 64 , 

with Sn^ = nj 1 — 8, then correspondingly 



where 



(S c denotes the area of a loop C and a is the lattice constant) and 



(2.41) 



(2.42) 



(2.43) 



(2.44) 



(2.45) 



lec 



G 



So the dynamics of is determined by the fluctuations of the density of holons on lattice. But since spinons and 
holons here are treated as independent degrees of freedom, one may neglect the dynamical effect of SA^ on spinon part 
at the mean-field level and replace it by some random flux fluctuations with a strength per plaquette equal to 5<j) (One 
may estimate 6<p « 7ry (Sn h ) 2 ). This can be justified at low temperature when a Bose condensation of holons (which 
corresponds to a superconducting condensation as shown later) occurs, where SA^ is expected to be substantially 
suppressed. On the other hand, however, in the high temperature phase where the motion of holons is much less 
coherent, the fluctuation effect of <5A' l „- can dominate over A% and the separation of the latter from A' 1 , then becomes 
meaningless. In this case, one may approximately describe the effect of A^ as a collection of randomly distributed 7r 
flux quanta with the number equal to that of holons. In both limits, the dynamics of A\: may be neglected. 



Including the hopping term Ht 



First of all, we note that the wavefuntion w ma (i) as the solution of the linear equation (2.32) is not unique, and it 
can be always multiplied by an arbitrary global phase factor e laXm : i.e., 

^(i) - e^-«w(i), (2.46) 

without changing the order parameter A|- and the mean-field state. Correspondingly the Bogoliubov transformation 
can be generally rewritten as 

ha = ^2 \ U mlma - ^mlm-a) s""* W ma (l) • (2.47) 



In particular, e lCTXm can depend on the holon configurations because the Hilbert space of bi a is only well-defined at 
each given holon configuration due to the no-double-occupancy constraint. The hopping term will mix the Hilbert 
space of bi a at different holon configurations, and such a freedom in phase choice can be fixed by optimizing the 
hopping integral of holons below. 



Now consider the hopping term H t in (2.5). By using the Bogoliubov expression for spinon operators, a straight- 
forward calculation gives 



e^vbUu, ) = V e ^< CT 0> mCT (i) e -- A ^ [vl + (u 2 m + !&)< 7 L7™>] • (2-48) 



\ a I ma 



Note that A^ m in the above expression denotes the difference of Xm before and after the holon changes the position. 
If one simply chooses A\ m = 0, namely, the phases of w m a{i) to be the same for all hole configurations, then 
the right-hand side of ( |2.4£ ) vanishes for the nearest-neighboring i and j. This can be verified by noting that 



Wfna{i) = (— f) 4 Wmcr is also a solution of (2.32) with an eigenvalue ^ m — — £ m and the cancellation in (2.48) stems 
from the fact that those quantities like v? m and only depend on E m in ( [2 .46 ) which is symmetric under £ m — > — £„ 



But such a cancellation is removable by a simple choice of the phase shift in e lCTXm at different holon configurations 
when each time a holon changes sublattice sites: 

- [sgn(U)} x e l " Xm , (2.49) 

or 

e- laAx ™ = -sgn(U)- (2-50) 



Then, one finds 



B ° ^ £ (e^jAr) = ^ £ B °m K + « + OklLrlm*)} , (2.51) 



where B° m = J^m e laA ^w* ma {j)w m a(i) [~sgn(^ m )} is given by 



B° m = - fl$n(£m)£<«) e l ^w* ma (j)w mr7 (t) 
= sgn(^ m )§f- Y,i w* ma {i)w ma {i) 

%7 > ■ (2-52) 
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In obtaining the second line above, we have used (2.32) with Af ■ = A s and J, 



|AV. 



Holons thus acquire a finite hopping integral without introducing any extra order parameter. The effective holon 
Hamiltonian is given by: 



H h = -thj^e^hjh, + H.c, 

(ij) 



which is derived from Ht with an effective hopping integral 

th = tB . 



(2.53) 



(2.54) 



It is important to note that such a finite kinetic energy (th ~ t) that each holon has gained on the present mean-field 
spin background cannot be similarly realized in the slave-fermion-Schwinger-boson scheme, exactly due to the hidden 
phase-string effect: the sign a — ±1 in (2.3) will always lead to th = and make a spiral twist (a new order parameter 
in favor of hole hopping) necessary in any local mean-field treatment. 

Finally, to be consistent, the hopping effect on spinons is obtained from Ht as: 



Hi 



H.c. + AJ h B N, 



(2.55) 



in which Jh = (e lA ^h\hj)t oc St measures the strength of hopping effect on the spinon part. (The constant in ( 2.55 ) 
is introduced such that (Hi) — 0.) Hence the total Hamiltonian describing spinon degrees of freedom is composed of 
two terms 



H s — H' 



Hi, 



(2.56) 



where H J in ( 2.19| ) has been diagonalized at the mean- field level before. iJ* can be expressed in terms of (2.47) by 
noting that b^ a and bi„ in ( [2.55] ) should differ by a phase shift ( 2.49| ) as a holon switches sublattice sites. It then gives 



H l - -Jh^^Bmi 



r ) (^mT" 



H.c. + 4J h B N, 



(2.57) 



after using the orthogonal condition w^ a (i)w m i a (i) — 8 m ,m' as well as (2.32). On the other hand, one has 
H J = Ylmo- ErnJm a "fmcr+ const, (note that X m = X inside E m here). Then H s can be diagonalized in a straightforward 
way. If we are still to use E m to denote th spinon spectrum for the sake of compactness, then H s can be finally written 
as 



H x 



E 



(2.58) 



where 



But here the spinon spectrum 



£ S = -2j2 E mn(Er, 



J 



{ij) 



(2.59) 



(2.60) 



is different from the previous one obtained in previous section by a correction to A m due to the hopping effect: 

Jh 



A, 



A 



(2.61) 



Therefore, the hopping effect on the spinon part is solely represented by a shift from A to A m in ( p.61| ). The 
Bogoliubov transformation (2.61) remains unchanged and so do u m and v m defined in ( 2.35|) a nd (2.36), so long as 
the renormalized A m is used. The Lagrangian mutiplier A in (2.61) is still determined by (2.4C) where both A m and 
E m should be replaced by the renormalized ones in (2.60) and (2.61). Finally, the self-consistent equation ( 2.38Q for 
the RVB pairing order parameter maybe rewritten as 
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A* = i_M y J|_ C oth (2.62) 
AN ^ J s E m 2 v ; 

(Note that if the Bose condensation occurs, one may separate the contribution from the condensation part on the 



right-hand-side by A^ c = (1 — 25)\£o\ 2 n b gc /4:J s \o). In obtaining (2.62), we have used an approximate relation 
l/(2N)J2(ij) (( A W 2 ) h ~ A s2 /(1 - 2(5) at 5 < 0.5 limit. Such a relation can be obtained by assuming A?- = Af 
when i and j belongs to occupied sites and Afj = if i or j is at hole site and by noting that each hole accounts for 
Af,- = at four adjacent bonds at dilute-hole limit which leads to A s = (1 — 25) Af . 



III. PHASE DIAGRAM 

A. Unified bosonic RVB phase 

Our mean-field theory has been constructed based on a single bosonic RVB order parameter A s . Such an order 
parameter controls short-range spin-spin correlations in both undoped and doped regime. Fig. 3 shows a typical region 
of A s ^ obtained by solving the mean-field equations which has been briefly discussed in Ref. |Is|| . It obviously 
covers the whole experimentally interested temperature (from T = to T ~ 0.5 — 0.9 J/fcs) and doping (from 8 = to 
S > 0.3) regime. Several low-temperature regions within this phase as marked in Fig. 3, including the superconducting 
phase, will be discussed in the following sections. The normal-state within this phase will correspond to a "strange 
metal" phase, where magnetic and transport properties are expected to be different from conventional metals. It 
is noted that in the bosonic RVB description of spin degrees of freedom, the order parameter A s does not directly 
correspond to an energy gap, in contrast to the fermionic RVB theory JlJ (the latter is similar to the BCS theory in 
mathematical structure). Also note that the crossover from A s ^ phase to A s = phase at high temperature is 
similar to the half-filling case fill] which does not correspond to a real phase transition. 



In obtaining A s in Fig. 3 by solving (2.62), one also needs to determine the spectrum £ m from (2.32) and decide 
the chemical potential A in terms of (2.4C). We have chosen the parameter Jh — SJ (which corresponds to t ~ J) 
and solved £ m under = A^ , but other choices of Jh as well as including the fluctuating part SA^ do not change 
significantly the range covered by A s ^ 0. The effect of SA^ will be the subject of discussion in the next section, and 
we will always use the same Jh below. 

Although A s 7^ practically covers the whole doping regime, at a larger doping concentration, this mean-field theory 
may no longer be energetically favorable due to the competition between the hopping and superexchange energies. 
Actually, the phase string effect itself is an indication that the bosonic description of spins leads to frustration of the 
motion of doped holes, and vise versa. With the increase of doping concentration, one possibility is that eventually a 
statistical-transmutation may occur to effectively turn bosonic spinons into fermionic ones as to be discussed in Sec. 
IV. Beyond such a point, the present mean-field theory will break down, which may determine a crossover to the 
so-called overdoped regime. We will explore this issue elsewhere. 

B. Bose condensation of spinons: AF ordering vs. phase separation 

1. AFLRO and insulation phase 
At half-filling, the spinon spectrum E m is known to be gapless at zero doping and zero temperature which ensures 



a Bose condensation 20 of spinons. Such a Bose condensation of spinons, as represented by n h BC ^ in (2.40) 



describes a long-range AF spin ordering. [|20| The Bose condensation or long-range AF order can be sustained up 
to a finite temperature Tjv > if the three-dimensional effect (interlayer coupling) is included. In the following, we 
consider how this AFLRO picture evolves at finite doping. 

Based on the expression (2.13), spin operators, Sf and Si , can be easily written down in terms of spinon operator 
bi a after using the constraint (2.7): 



s* = ~$>4a- (3-1) 



and 
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S+ = ^ M _i)V** (3.2) 

and Sr = (fit)t. 

At 5 = 0, one has = 0, and the Bose condensation leads to 

(St) « ("!)*, (3-3) 

i.e., an AFLRO. But at 5^0, even when the spinons are Bose-condensed, (S^~) should generally vanish due to the 
fact that 

(e^*)=0. (3.4) 

The proof here is straightforward. Note that in the definition of <I>^ in ( 2.16| ), the angle ^(/) can be transformed 



6i(l) -» 6i(l) + </> (3.5) 
for an arbitrary cf> without changing A^, A{j, and thus the Hamiltonian. But e t$i [( 2.16| )1 changes accordingly 



x e^ N , (3.6) 



Here N h is the total holon number. Thus the average of such a phase must vanish at finite doping as given in 
Since <j>^ describes vortices centered at holons, it is like a free- vortex phase as holons move around freely in metallic 
phase, which resembles a disordered phase in a Kosterlitz-Thouless-type transition. 

Only in the case that holons are localized (i.e., in insulating phase) such that the phase string effect is ineffective, 
the AFLRO may be recovered. In this insulating phase, holons are perceived by spinons as localized vortices like 
in the mixed state of a type-II superconductor, and by forming "supercurrents" to screen those vortices, $^ in St 
can be effectively canceled out by the opposite vorticities generated from spinons. After all, the phase string effect 
is no longer important if the motion of holons is limited. We do not know if the localization of holons can happens 
intrinsically or has to be under some external factors like impurity effect. But we expect such an insulating phase to 
exist only at a very dilute density of holons at the expense of latter 's kinetic energy. 



2. Bose condensation of spinons in metallic phase: underdoping 

We have shown that the AFLRO must be absent in metallic phase. One may naturally wonder if the Bose conden- 
sation of spinons can still persist into metallic phase, and if it does, then what is its physical meaning? 

To answer these questions, let us first to inspect how the spinon spectrum E m is modified by doping. In Fig. 4, we 
compare the spinon density of states p s (E) at 6 = 1/7 « 0.143 (solid curve) with the 5 — case (diamond curve in 
the insert). Here p s {E) is defined by 

Ps(E) = ^Y.^ E - E ^)- ( 3 - 7 ) 

TCI 

One notices that a unique peak-structure is clearly exhibited at S = 0.143. This can be easily understood by noting 
that the spinon spectrum E m is basically determined by £ m which, as the solution of ( 2.40| ), has a Hofstadter structure 



(or the Landau levels in the continuum limit) due to a uniform flux (0 = nS per plaquette) represented by the vector 
potential threading through the square lattice. The broadening of the solid curve in Fig. 4 is due to the 
redistribution of eigenstates under the fluctuating flux SA^ , which is treated as a random flux (in white noise limit) 
here with a maximum strength chosen at Sep — 0.30 per plaquette. By contrast, the dashed curve marks the positions 
of sharp peaks of density of states in the limit of Sep = 0. 

We find that the Bose condensation of spinons can still occur at 6 = 0.143 with 8(p = 0.30 (but not at Sep = 0) 
as the solution of (2.4C). Recall that the spinon Bose condensation stems from the equation ( p. 40 ) with a nonzero 



n b BC representing the density of spinons staying at E m = state. In this case, there would be no solution at low 
temperature unless A takes a value to make E m gapless such that n b BC ^ can balance the difference between the 
left and right side of the equation, similar to the half-filled case. |2(J In particular, such a Bose condensation is found 
to be sustained up to a finite temperature Tbc ~ 0.21J even in the present 2D case. This is due to the vanishingly 
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small weight near E = in the density of states (in Fig. 4, there is a small tail in the solid curve which extends 
to E = 0), where the spinon excitations at low-temperature are not sufficient to destroy the Bose condensation as 
first pointed out in Ref. pi] ]. It is noted that in principle the strength 5<f) of the random fluctuations of SA^ should 
be self-consistently determined by the density fluctuations of holons. But here for simplicity we just treat 5<p as a 
parameter and then study the qualitative characteristics under different values of 5(f). The actual strength of the 
fluctuations of will be only crucial in determining the location of the phase boundary. 

Phase separation Note that the Bose condensation means a thermodynamic number of spinons staying at E m — 
— the lowest energy state which corresponds to the band edges of the spectrum £ m . So due to the Bose condensation, 
such quantum states will acquire a macroscopic meaning. But these band edge states of £ m are very sensitive to the 
fluctuations of A^ and the density of holons. Physically, the density of holons is fluctuating in real space (which is the 
reason leading to the fluctuations of A^) and there always exist those configurations in which the density of holes are 
relatively dilute in some areas where the flux described by A^ is reduced such that the band edge energies of £ m can 
be close to ±2 J s . Since the probability would be small for such kind of inhomogeneous hole configurations, the density 
of states generally looks like a tail — i.e., the Lifshitz tail — near the band edges. Hence, the corresponding Bose- 
condensed state will generally look like having a charge inhomogeneity, or, phase separation, with spinons condensing 
into hole-deficient regions to form short-range spin ordering. The true AFLRO is absent here. 

Pseudo-gap behavior The Bose condensation of spinons will also lead to a pseudo-gap phenomenon. In magnetic 
aspect, for example, the density of states shown in Fig. 4 indicates the suppression of spinon density of states between 
zero and the lowest peak, which stabilizes the Bose condensation as mentioned earlier. Since in the Bose condensation 
case there must be some residual density extending to E m = 0, a pseudo spin gap is thus present in p s (E). Its effect 
in the dynamic spin susceptibility will be discussed later. In the transport aspect, holons which are the charge carriers 



are scattered off by the gauge field A? A according to the holon effective Hamiltonian (2.53). Anomalies in transport 
properties have been found in Ref. [pl| with interesting experimental features in the similar effective Hamiltonian, 
where fluctuating fluxes depicted by A\^ play a central role. But the Bose condensation of spinons will lead to 
a substantial suppression of A|- and thus a reduction of scattering to holons. Hence, the Bose condensation also 
provides an explanation for the so-called pseudo-gap phenomenon shown in the underdoped high-T c cuprates, where 
the transport properties deviate from the high temperature ones below some characteristic temperature scale. 

Therefore, if the Bose condensation of spinons happens in metallic phase, it will result in a phase-separation or 



spin pseudo-gap phase without the AFLRO. With the increase of S, the reduction of the left hand side of (2.40) will 
eventually make the Bose condensation term, if exits, disappear from the right hand side. So the Bose condensation 
in general may only exist at small doping regime, which can be defined as the "underdoping" regime. In Fig. 3, the 
shaded curve sketches such a region outside the true AFLRO phase which is in a much narrower region (dotted curve) 
at finite doping. 

C. Superconductivity 

In the phase-string representation, the operator of superconducting order parameter 

Ag C = $>c*Ci-* ( 3 - 8 ) 



can be expressed in terms of ( 2.13 ) as follows 



in which 



(ftJeW) (*}e*«5) (-I)*, (3.9) 



a;, Y^, - rV h„-h, ,. (3.io) 

This is the basic expression to be used in the following discussion of superconducting condensation. 



1. Mechanism 



For the sake of simplicity, we will focus on the nearest-neighboring pairing with i — nn(j) below. Since the whole 
mean-field phase is built on 



11 



(A£)=AVO, (3.11) 
the electron pairing order parameter g AfP = (AfP) can be written as 



((ftJeW)(fct c l«5)) ( _i)i. (3.12) 



We see that the spinons are always paired in the present phase, as described by A s , up to a temperature scale ~ J at 
small doping. Thus, in order to have a real superconducting condensation below a transition temperature T c , the holon 
part has to undergo a Bose condensation or, strictly speaking in 2D, a superfluid transition in the Kosterlitz-Thouless 
sense (recall that both spinon and holon are bosonic in the present representation). 

One may notice that this superconducting condensation picture is somewhat similar to that in the slave-boson 
mean-filed theory. Q But there are two crucial differences. 

Firstly, the spinon pairing in the present case practically covers the whole superconducting and normal-state regime 
that we are interested in. In other words, A s ^ in the present mean-field theory defines a "strange" metal, and 
the normal-state anomalies ol experimental measurements in the cuprates, including the magnetic properties and 
the transport properties, are all supposed to happen within such a phase. In contrast, in the slave-boson mean-field 
approach the fcrmionic spinon pairing is directly related to the gap in the spinon spectrum and has to disappear at a 
much lower temperature scale beyond which "strange" metallic properties presumably start to show up. 

Secondly, the Kosterlitz-Thouless transition temperature of holons are believed to be much higher than the real 
T c in the cuprates, and thus one has to introduce other mechanism (e.g., gauge-field fluctuations pjj) to bring down 
the temperature scale in the slave-boson mean-field approximation. On the other hand, there is a unique feature in 
( |3.12 ), namely, the presence of phases like e^* ; . Here represents a structure of vortices (anti- vortices) centered 
at | (|) spinons. At T = 0, all f and J. spinons are paired up at a finite length scale and so are the vortices and 
anti-vortices in <&\, which implies Afp ^- as long as holons are Bose-condensed. At finite temperature, even though 
the Kosterlitz-Thouless transition temperature for hard-core bosons can be much higher, the "phase coherence" in 



can be more quickly destroyed at a lower temperature due to the dissolution of the vortices and antivortice 



bindings in <!>J after free spinons appear. Here the argument for (e^^e 5 *^) = is similar to the previous one for 
(e 2 ** ) = which corresponds the disappearance of the long-range AF order once holons become mobile in the metallic 
phase. This provides us an estimate of the upper limit for the superconducting transition temperature T c below. 



2. An estimate ofT c 



The holon effective Hamiltonian Hh in (2.53) determines the interaction between holons and those vortices described 
by A^j. If free vortices are few, the condensed holons may easily "screen" them by forming supercurrent, which will 
then effectively keep AfP finite. But if the number of free vortices, or excited spinons, becomes comparable to the 
number of holons themselves, one expects that the "screening" effect collapses and thus Afp = 0. It predicts that T c 
will be basically determined by the spinon energy scale in the following way: 



2 £ ^n(E„ 



kNS, 



T=T C 



1 and the left-hand side represent the average number of excited spinons in Ylia ^ia^is. 



(3.13) 



with the Bose- 



where n 

condensed part (if exists) excluded. The dashed line in Fig. 3 represents the T c 's determined by ( 3.13 ) in the limit 
S(f> — > 0. This curve may be regarded as represents the "optimized" T c , because with introducing the flux fluctuations, 
56 0, there is a finite density of states of spinons emerging at lower energy which effectively reduces T c defined in 
(3.13). In the "optimized" limit of 8<p = 0, one may further simplify Eq.(3.13) by only retaining the contribution from 
the lowest-peak (which has a degeneracy 6N/2) and obtains: 

-E„ (3.14) 



where c is given by 



In 1 + -V1 + W > 1 



(3.15) 



Here E s and £ s are the energies of E m and £ m , respectively, corresponding to the lowest-energy peak shown in Fig. 4 
at S(j) = 0. Therefore, T c is indeed determined by the characteristic energy E 8 of spinon excitations. 
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3. d-wave symmetry of the order parameter: phase string effect 



Finally, let us briefly discuss the symmetry of the order parameter A^p. Basically, one needs to compare the relative 



phase of A?- c between j = i + x and j = i + y, or the phase change of the quantity h/. = /ije 2 *j in ( 3.12 ). Imagine 
that we move a holon from j = i + x to i + y via site i + x + y. At each step the holon has to exchange positions with a 
spinon with index ay at site j 1 which leads to an extra phase ay = ±1 due to e%®i in h .. Even though other spinons 
outside the path also contribute to, say, <& b i+x+y — 3>£_|_£, but their effect is canceled out as hj picks up the same phase 

change but with opposite sign in H^. Therefore, in the end hj acquires a total phase Oi+y ■ Oi+x+y which is just the 
phase string on such a path. Its contribution is always negative on average for a short-ranged AF state. Assuming 
that this is the dominant path, one then concludes that A^p has to change sign from j = i + x to i + y, namely, 
the d-wave symmetry. If only the nearest-neighbor-site electron pairing is considered, the order parameter in the 
momentum space can be written in the form A sc (k) cx (cos k x a — cos k y a). Therefore, the nonrepairable phase string 
effect and AF correlations are directly responsible for the d-wave symmetry of the superconducting condensation in 
the background of A s ^ 0. 



D. Experimental implications: dynamic spin susceptibility 

1. Local spin dynamic susceptibility 

The dynamic spin susceptibility function Xl( w ) = V-^ Xzz(*; *> w ) describes the on-site spin dynamics and is 
derived in Appendix as follows 



771771 ' \ i(7 



mm' > 

+ n (E m ) + n(E' m ))(u m v m , + v 2 m u 2 m ,)6(\u>\ - E m - E m ,) 

+ (n(E m ) - n(E m ,))(u 2 m u 2 m , + v 2 m v 2 m ,)8{u + E m - E m ,)] , (3.16) 

where the summation only runs over those m's with £ m < (note that E m is symmetric under £ m — > — ^ m ). If 
there is a Bose condensation of spinons, the contribution from the condensed part to x'l m ay be explicitly sorted out 
as 



X, 



(u) = sgn(uj) (l^sc) ^ E K 0m^ !L S(u) - E m ), (3.17) 



with Ko m = N Y)j„ \w o a (i)\ 2 \w ma (i)\ 2 where the subscript refers to the E m — state. 



Based on Eq.( 3.16 ), two kinds of mean-field solutions, with and without Bose condensation of spinons, will be 
studied below. Without loss of generality, we consider these two cases at S = 0.143 in which the corresponding density 
of states, p s (ui), is already shown in Fig. 4 for both cases. Let us firstly focus on the lowest peak of p s (E) shown in 
Fig. 4. The contribution of such a peak to Xl( w ) was previously discussed in Ref. [is] and is illustrated in Fig. 5 by 
the lowest sharp peak (dashed line) for 8<fr — and the lowest twin peaks (solid curve) for 5(f> = 0.30, respectively. 
One sees two very distinct features here. For the case of 5<f> — 0.3</>, there is a Bose condensation contribution at 
T < Tbc ~ 0.21J and it leads to a double-peak structure. But in 5<f> = case, the Bose condensation is absent and 
one finds only a single sharp peak at 2E S ~ 0.4J. In other words, Xl( w ) has drastically different characteristics for 
cases with and without a spinon Bose condensation. 

The twin-peak splitting may be understood as follows. The second peak corresponds a pair of spinons excited from 
the RVB vacuum, while the first peak describes a single spinon excitation as the other branch of spinon is in the 
Bose-condensate state. Such a lowest-peak basically maps out the lowest peak of the spinon density of states p s (E) in 



Fig. 4 according to (3.17), and the second peak in x'l( w ) is located at an energy approximately twice larger than the 
first one. The latter one will be always around at low temperature no matter whether there is a Bose condensation 
or not. So there is a distinct behavior of those two peaks at different temperature as shown in Ref. fTsfl , where the 
weight of the lowest one gradually diminishes as the temperature approaches to Tbc- 

In contrast, there is only one sharp "resonance-like" peak left in the where the Bose condensation is absent. It 
corresponds to a pair of spinon excitations located at the lowest peak (E 8 ~ 0.2 J) of p s in Fig. 4 with 5(f) = (dashed 
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curve). Note that this is the limit where the flux fluctuating part 5 is totally suppressed such that a real spinon 
gap is opened up at low-energy as shown in Fig. 4. The location of the "resonance" peak, 2E S , is slightly lowered in 
energy as compared to the corresponding (second) peak in the twin-peak case at 5(f) — 0.3</>. The energy scale of this 
peak (2E S ~ 0.4J) at 5<fi — > limit is roughly independent of Jh and thus of t. This is because Jh term only shifts A 
to A m by a constant if there is no dispersion in £ m near the lowest (highest) peak, which will not affect E m near the 
lowest peak as A will readjust its value correspondingly. 

2. Underdoping vs. optimal- doping 

We have previously discussed the Bose condensation of spinons and argued that it exists only in an "underdoped" 
regime. We have also shown that holons as bosons can experience a Bose condensation at T c , leading to the su- 
perconducting condensation. If the Bose condensation temperature for holons is higher than Tbc for spinons, i.e., 
T c > Tbci the holons will become Bose-condensed before spinons do and it will be generally a uniform state since 
there is no inhomogeneous spin ordering above Tbc- Thus one may have w A^ with Sip being much less than in 
the normal state. In this case, the Bose condensation of spinons can be effectively prevented at lower temperature 
because a homogeneous A^ generally leads to an opening of a real spinon gap as shown in Fig. 4 for <j> = 0. So to 
be self-consistent, once T c > Tbci Tbc may no longer exist and T c becomes the only meaningful temperature scale. 
Furthermore, we have already seen that T c is also optimized under 6cj> ~ 0. Therefore, this region may be properly 
defined as the "optimal-doping" regime in our theory in contrast to the previously defined underdoped regime at 
Tbc > T c . 

Such an optimal-doping phase is charge homogeneous and the Bose condensation of spinons is absent. It is charac- 
terized, below T c , by the "resonance-like" peak emerging in x'Li 1 ^) & t m Fig- 5- It is in accord with the AlmeV 
peak found |22| in the optimally doped Y BaiCu-zOi below T c , if one chooses J = lOOmeV here. Above T c , the 
"resonance" peak will quickly disappear as the motion of holons becomes incoherent and a different behavior of A^ 
is involved as discussed in Sec. II. 

On the other hand, the underdoping regime with Tbc > T c is characterized by a low-energy twin-peak structure 
in at T < Tbc- In contrast to the "optimal-doping" case, such an energy structure may not be qualitatively 

changed even when T is below T c , as holons are also expected to be condensed inhomogeneously in favor of the spinon 
energy. A twin-peak feature has been observed recently in the underdoped YBazCu^Os^ compound by neutron 
scattering [ p3| in the odd symmetry channel. In the experiment, the lowest peak is located near 30meV and the 
second one is near 60meV, indeed about twice bigger in energy. Most recently, in the underdoped Y Ba2Cu^OQ.Q, a 
second energy scale near ~ 70meV has been also indicated besides the earlier report of the lower energy peak near 
34:ineV. p5[ ] It is noted that the energy scales shown in Fig. 5 are generally doping dependent, and those energies in 
Y BaiCu-zOss are expected to be relatively smaller than the corresponding peaks in Y BaiCu?,0§.§. 

A word of caution about the comparison with the YBCO compound is that the latter is a double-layer system 
where two adjacent layer coupling is also important. But we do not expect the double-layer coupling to qualitatively 
change the above energy structure of x'l i n the odd symmetry channel. We point out that the fluctuating part of the 
gauge- field SA^ usually makes the two adjacent layers difficult to couple together unless there are AF spin domains 
in the charge deficient region where the total A^ is suppressed, as may be the case of phase separation. But in the 
uniform phase, with 6A^ being suppressed below T c due to the Bose condensation of holons, the effective coupling 
between layers can also be greatly enhanced to gain interlayer-coupling spin energy. The Anderson's confincmcnt- 
deconfinement phenomenon p6| may become most prominently in the optimal-doping regime, which needs to be 
further explored. 

3. Prediction 

Fig. 5 also shows Xl( w ) i n the whole energy regime at 8(f>/4> = (dashed curve) and 0.3 (solid curve), respectively. 
As compared to the 6 — curve in the insert, a multi-peak structure is present as well at high energies for such a 
doped case. For example, if the AlmeV peak in YBa 2 Cu 3 07 is explained by the lowest peak in the case of 8<j) = 0, 
then the theory predicts a second "resonance" peak near 120meV to be found. Those high energy peaks in Fig. 5 
become rather closer in energy especially in the Bose condensed case (solid curve) which could be very easily smeared 
out either by the experimental solution (it may be further complicated by the fact that the momentum-dependence 
varies drastically among those peaks) or by the dynamic broadening due to the finite life time of spinons which is 
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beyond the present mean-field treatment. Nevertheless, the multi-peak structure, especially the twin-peak feature 
at low energy in the spinon Bose-condensed case, should become observable by high solution measurement at low 
temperature as the unique prediction of the present theory. 



IV. DISCUSSIONS 



In this paper, we have approached the doped antifcrromagnct from the half-filling side, where the bosonic RVB 
description is known to be very accurate for the antiferromagnetism. The crucial modification at finite doping comes 
from the phase string effect induced by doped holes. For example, doped holes are turned into bosonic holons by such 
a phase string effect so that both the elementary spin and charge excitations are bosonic. The Bose condensation 
of spinons in insulating phase and the Bose condensation of holons in metallic phase determine the AFLRO and 
superconducting phase transitions, respectively. 

While the bosonic RVB pairing, representing short-range AF correlations, is always present and is the driving 
force behind the antiferromagnetism and superconductivity, it is the combination with the phase string effect that 
decides when and where they occur in the phase diagram. For instance, the Bose-condensation of spinons leads to the 
AFLRO only in the case that holes are localized. In the metallic phase where holes become mobile, the AFLRO will be 
destroyed by the phase string effect. But the Bose-condensation of spinons may still persist into weakly doped metallic 
region, leading to an "underdoping" metallic phase with charge inhomogeneity (phase separation) and pseudo-gap 
phenomenon. 

There still are many theoretical and experimental issues which have not been dealt with in the present paper and 
are left for further investigation. Here we conclude by giving several critical remarks. The first is about the phase 
diagram at larger doping. Recall that in our mean-field description, the metallic phases are characterized by two 
temperature scales: T c and Tbc, and we have argued that T c > Tbc determines an "optimal-doping" regime at low 
temperature where Tbc is no longer meaningful. But beyond this regime, there is a possibility that holons may tend 
to be always Bose-condensed even at normal state in favor of the hopping energy. If this occurs, the gauge field in 
Hj may have to be "expelled" to the spinon part, leading to a statistics transmutation to turn spinons into fermions 
and causing a collapse of the bosonic RVB order parameter at the normal state. In this picture, the normal state in 
overdoped regime may simply recover the fermionic uniform RVB state. Q] 

The second issue is about the time-reversal symmetry. Recall that the sharp peak structure in spinon spectrum 
below T c is mathematically similar to a Landau level structure in a uniform magnetic field. One may naturally wonder 
if some kind of time-reversal symmetry would be apparently broken like in the anyon theories. || Below we point out 
that this is not the case in the present theory. First of all, it is easy to see that there is no breaking of the time-reversal 
symmetry in the holon Hamiltonian Hh (2.53) in which the gauge field AL = Afj — <j>1j. Here Afj behaves like a 
fluctuating gauge field with (Afj) = 0, and <£°- describes a uniform 7r-flux per plaquette which does not break the 
time-reversal symmetry either as a gauge transformation can easily change n flux into —n flux per plaquette. As for 
the spin part, even though spinons see A^ which breaks the time-reversal symmetry, one should remember that the 
physical observable quantity is the spin-spin correlation functions like x'L shown in ( 3.16| ), which can be easily shown 
to be invariant under A^ — > — A\, . 

Lastly, the sharp peak in the spinon spectrum in the uniform phase provides an explanation for the AlmeV peak in 
the neutron-scattering measurement of YBaCuO "90 K" sample, but it also means a real gap in the spinon spectrum. 
From a naive spin-charge separation picture, one would expect the single-electron Green's function to be a convolution 
of spinon and holon propagators and the electron spectrum may also show a finite gap as well in the superconducting 
state which may be inconsistent with d-wave symmetry. But we note that in the present spin-charge separation 
formalism ( 2.13 ), there is an additional phase field Q s ^ n9 representing the phase string effect. It means that if a 
"bare" hole created by Ci a decays into mobile spinon and holon, a nonlocal topological effect will be left behind which 
could cost a logarithmic-divergent energy. In other words, the phase string effect in 2D case will serve as a confinement 
force to prevent a newly doped hole from dissolving into elementary excitations. (Of course, internal charge and spin 
excitations without involving the change of total electron number are still described by the spin-charge separation in 
the present mean-field theory.) In fact, even in ID case, the single-electron Green's function looks quite differently 
from a simple convolution of spinon and holon propagators, and recently Suzuura and Nagaosa |2?j have discussed 
the crucial role of the phase string effect in understanding the angle-resolved photoemission spectroscopy in SrCuC>2 ■ 
p8[ Our preliminary investigation in 2D indicates that a bare "hole" wavepacket injected into the background of 
the spin-charge separation mean-field state will behave more like a conventional band-structure quasiparticle in a 
Fermi liquid which shows d-wave gap structure when the holons are Bose-condensed and pseudo-gap structure when 
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spinons are Bose-condensed. So experiments involving injecting electron or hole into the system like photoemission 
spectroscopy may no longer provide direct information of elementary excitations like in the conventional Fermi liquid 
theory, due to the confinement of phase string effect. 
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APPENDIX A: DYNAMIC SPIN SUSCEPTIBILITY FUNCTION 



Local spin susceptibility function is defined in the Matsubara representation as follows 



-<T T Sf(r)Sf(0)>, 



(Al) 



where io n — ^jp-- In the following we will determine the dynamic spin susceptibility function x'Lpi^ *! w ) based on the 
present mean-field theory. 

Consider (T T Sf (t )Sf (0)) in a = (3 — z case. In the present mean-field formulation, one has 



(A2) 



(T t S?(t)S?(0)) = ^aa'(T T 6i ff (r)6 iCT (r)6t(0)6 iCT ,(0) 
= \ E (2V&t»M0)) (TrhAr)bl(O) 



In terms of ( 2.47 ) and ( |2.5S| ), one has 



(A3) 



Then by noting w* m _ a = w ma and that for each m with £ m < 0, one always can find a state m with = — £ m > 
with a wavcfunction 



Wfhcr{i) = {-l) % w mcr (i) 



(A4) 



according to (2.32), we get 

(T t C(t)MO))t>o = w* ma (i)w m<J (i) x [uln(E m )e E ™ T +vl(l + n(E m ))e- E ™ T ] , (A5) 

m 

where the summat ion ^2' m oru y run s over those states with £ m < 0. Then it is straightforward to obtain Xzz after 
integrate out r in (Ia5): 



(A6) 



where Xzz is defined by 



X ( z f(i,i]iu n ) = ^E K w(*'*) 



/ ± ^n{E m ,) - n(E, 



+ (ti m ') 2 (l + n(E m )+n(E m ,))^ 



1 / 1 



iuj n + E m — E, 
1 



2 V iujn + E m + E m i iui n — E m — E m i 



(A7) 
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with 

= J2 \ w ^)\ 2 l^mv(i)| 2 ■ (A8) 



Here the coherent factors, £> mm / and l mrn t are denned by 



mm 
rum 



U"mV m ' lb V m ll m ' . yA.Q) 



Finally, the dynamic spin susceptibility function y^' zz (i^i\Lo) can be obtained as the imaginary part of after an 
analytic continuation ioj n — > u; + is made: 

i; w) = i; w) + *<+>(*, i; w), (A10) 

where 

mm' 

•£(M - £ m - £ m ,) + 2 (n(£ m ) - n(E m 0) (Pw) 2 <*(w + E m~ E m ,)] . (All) 
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Fig. 1 A sequence of sign mismatches (with reference to a spin background satisfying the Marshall sign rule) left 
by the hopping of a hole from site a to site b on square lattice. 



Fig. 2 Fictitious n flux-tubes bound to holons which can only be seen by spinons and are described by the gauge 
field A% defined by (pig). 



Fig. 3 Phase diagram of the bosonic RVB state characterized by the order parameter A 3 (solid curve). Within 
this phase, the shaded curve sketches a region where a Bose condensation of spinons (BC) occurs, which leads to 
the AFLRO in an insulating phase (dotted curve) but a charge inhomogeneity with short-ranged spin ordering in 
metallic region which defines an underdoped regime. Superconducting condensation (SC) happens due to the Bose 
condensation of holons and T c (dashed curve) is determined under an optimal condition (see the text). 

Fig. 4 Spinon density of states at 6 = 0.143 and T = 0. Solid curve corresponds to 8(f) = 0.3(f) and the dashed curve 
is for 8(f) = 0. Insert: the density of states at half-filling. The energy is in units of J. 

Fig. 5 Local dynamic spin susceptibility versus energy at 8 = 0.143. Solid curve: 8(f) — 0.3(f) and dashed curve: 
8(f) = at T = 0. Note that the lowest peak of dashed curve splits into a twin-peak structure in solid curve (see text). 
The half-filled case is shown in the insert for comparison. 
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